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SUMMARY 


A numerical procedure has Been developed for analyzing wings and wing- 
Body comBinations and for designing optimum wing camBer surfaces in the 
presence of a Body. The method is very general and applies to wings of 
arbitrary planform and Bodies of arbitrary cross section and camBer. The 
P f° c f^ u f e has been Programed for automatic computation and considerable 

made to allow the user to analyze a great variety of con- 
figurations with relatively simple input data. 

For a given wing or wing-body combination, five classes of problems 
may Be solved: 


(1) Wing warp required to support a given loading on the wing 

(2) Wing loading for a given wing warp 

(3) Pressures on the upper and lower surfaces of a warped wing of 

small But finite thickness 

(4) Minimum drag wing shape for a given lift constraint 

(5) Minimum drag wing shape for a given lift and moment constraint 



INTRODUCTION 


The work described in this paper has Been undertaken in an effort to 
develop a procedure for designing wing-body comBinations with low drag at 


Several methods have been published which enable one to compute the 
camBer surface of minimum drag for an isolated wing at a given lift. None 
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of these methods, however, considers the effect of the flow disturbance gen- 
erated by the body on the shape of the optimum camber surface. The new method 
used in this work enables one to determine the lifting surface of minimum drag 
for a given value of lift in the presence of a body which may be at incidence 
relative to the free stream. In addition to this computation, the lift and 
drag of a configuration of given geometry may be determined, thereby providing 
the wing-body interference for configurations of arbitrary planform and camber. 
By using this procedure, one can assess the effects of such modifications as 
wing and body camber and incidence on the lift, drag, and moment characteris- 
tics of a given wing-body configuration. 

The numerical procedure has been programed for automated computation for 
a great variety of configurations* 

The theoretical analysis and the computer program used for the numerical 
computations in this paper were developed by the Aerodynamic Research Unit of 
the Airplane Group of The Boeing Company under NASA contract NAS2-2282. (See 
refs. 1 and 2. ) 


SYMBOLS 


AR aspect ratio 

b span 

c chord 

Cj^ slope of lift curve, dCj^da, per radian 

C L lift coefficient 

C m pitching-moment coefficient 

Cj^ pitching moment at zero lift 

AHp difference in pressure coefficient between upper and lower surfaces 

of wing 

d body diameter 

L/D lift-drag ratio 

M Mach number 

P difference between pressure coefficient at given angle of attack and 

pressure coefficient at zero angle of attack 


410 



r 

body radius 

u„ 

free-streara velocity 

x,y,z 

Cartesian coordinates 

X* 

distance from the leading edge 

a 

angle of attack 

3 = {iF 

- 1 

A 

leading-edge sweep angle 


METHOD OP ANALYSIS 


The actual wing-body combination is replaced by distributions of singular- 
ities, which satisfy the linearized equation of supersonic flow, whose strengths 
are adjusted to satisfy the boundary conditions required by the geometry of the 
particular configuration. The wing warp and incidence are represented by dis- 
tributions of vorticity corresponding to pressure differences across the wing. 
The wing thickness effect is given by sources and sinks located over the wing 
reference plane. The body thickness, camber, and incidence are sim ula ted by 
line sources and doublets placed along the body reference axis. Finally, the 
interference effects of the wing on the body are cancelled by a distribution of 
vorticity on the surface of the body. 

A typical wing-body combination is represented in figure 1. A grid which 
conforms to the general geometry of the configuration subdivides the wing and 
body into a large number of small panels. Each panel is a region over which a 
particular singularity strength is held constant. In this way, the continuous 
distribution of singularities representing the actual wing and body is approxi- 
mated by simple functions with a finite range of values. Associated with each 
panel is a control point for matching the boundary conditions. At each control 
point, the flow must be tangential to the panel surface. In contrast to the 
wing-alone programs in common use which usually match only the downwash, the 
normal vector to the panel surface is used in this computation. The resultant 
normal velocity at each panel control point may now be expressed as a system of 
linear equations in terms of the singularity strengths. The coefficients of 
this system of equations are computed and stored as a matrix of aerodynamic 
influence coefficients. Then, for a given camber shape or pressure distribu- 
tion, the various singularity strengths which satisfy all the boundary condi- 
tions are computed, and from these, the pressure distribution, lift, drag, and 
moment on the wing and body may be computed. 

This method differs from that of Harris discussed in paper no. 27 in that 
the forces are computed by integration of the local pressures over the surface 
of the configuration. This method may be referred to as a near-field theory. 
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whereas that of Harris may he called a far-field theory. Of course, the results 
of the two methods will he identical except for the leading-edge suction force 
(on subsonic leading edges) which is not included in the surface- integration 
method. 

Since the lift and drag of the configuration are functions of a finite num- 
ber of real variables, the determination of the surface of min imum drag for a 
given lift constraint can now be solved by the Lagrange multiplier method. 

A complete description of the theory and programing details is contained 
in references 1 and 2. 


RESULTS AND DISCUSSION 


Pressures and Forces on Flat Wings 

The accuracy of the technique for computing pressures and forces on iso- 
lated wings is indicated in figures 2, 3, and 4. In each case, planforms have 
been selected for which exact solutions are known to the linearized equation 
of supersonic flow. 

In figure 2, the program results are compared with exact conical flow, 
theory for delta wings. Examples are shown for both subsonic and supersonic 
leading-edge wings. The symbols represent the pressures given by the program 
for wing chords located at 15, 45, and 75 percent of the semispan. The solid 
curves were computed from conical flow theory and may be found as solutions 3 
and 6 in reference 3* The results agree quite well except at the sharp ridges 
of the supersonic solution. This rounding of the pressures computed by this 
program occurs because the pressures shown are actually average pressures over 
a wing panel. 

The next example, illustrated in figure 3> a test of the ability of 
the method to predict tip effect and to account for a subsonic trailing edge. 
The planform is a constant chord wing with PAR = 1.92. The leading and 
trailing edges are subsonic with p cot A = 0.6. The program results along 
four different wing chords at 25, 50, 75 > and 95 percent of the semispan are 
given by the symbols. The solid lines are the results given in reference 4 
which were computed by the technique of superposing conical flows. As in the 
previous example, the program results agree very well except in the regions 
where the pressure distributions vary rapidly. Along the chords where the 
analytic solution varies rapidly, the program results underpredict the pres- 
sures in front of a sharp pressure drop and oveipredict those behind the drop 
with the result that the integrated loadings are quite close to those of the 
analytic solution. This can be seen by the results shown in figure 4. The 
integrated loadings for a delta wing with 60° of leading-edge sweep are shown 
by the symbols. The solid curve is the well-known exact conical flow solution 
for a delta wing. 
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Wing-Body Loadings 


To confirm the accuracy of the method for calculating wing-hody interfer- 
ence, a configuration whose loading has been measured was studied. A sketch of 
the configuration is presented as figure 5 • The rectangular wing of aspect 
ratio 3 is mounted in the midwing position on a cylindrical body with an ogival 
nose. The wind-tunnel model was constructed with pressure taps along 0° and 45° 
meridian lines on the body and along several chords on the wing. The pressures 
measured at a Mach number of 1.48 and reported in reference 5 are compared with 
the program results in figure 6 . The parameter P represents the value of 
pressure coefficient at angle of attack less the pressure coefficient at zero 
angle of attack. The predicted pressures agree quite well with those measured 
experimentally at a = 2°. As the angle of attack increased, the agreement 
was less satisfactory. 

The examples shown and other test cases studied provide confidence in the 
a ^^^-4ty of the method to predict loadings, forces, and moments to a degree of 
accuracy well within the requirements of engineering design. 


Prediction of the Effect of Wing and Body Camber 

All the preceding results are accessible by established methods and serve 
to provide confidence in the method. The results to be shown in this section 
are not within the scope of any of the well-known methods. 

The aerodynamic characteristics of a wing-body combination are sensitive 
to modifications, such as wing incidence relative to the body, wing or body 
camber, etc. In fact, it is precisely through modifications such as these 
that the characteristics of an airplane are brought in line with the design 
requirements. This computational method enables the designer to make accurate 
estimates of the effect of these changes to the basic configurations. 

As reported in reference 6 , a study was made of a series of wing-body wind- 
tunnel models employing various combinations of wing and body warps. The wing 
was arrow shaped with a leading-edge sweep angle of 70 ° and an aspect ratio of 
2.24. Configuration 1 consisted of a flat wing on the uncambered body. Con- 
figuration 2 was made up of a cambered wing mounted on the uncambered body such 
that the body axis was at zero incidence when the wing was at the design con- 
dition ) namely, C^ * 0.08 and M * 2 . Configuration 3 had the same cambered 

hut the axis of the body was alined with the root chord of the wing j also, 
the nose of the body was drooped and the rear portion of the body was swept 
upward. 

The effects of these variations in model geometry are illustrated in fig- 
ure 7 . The symbols represent the experimental data reported in reference 6 
and the curves are the program results. Since the theoretical method described 
in this paper does not include viscous effects, it was necessary to estimate 
the skin-friction drag in order to compare theoretical and experimental drag 
values. For this study, the skin friction was taken to be the difference 
between the experimental drag and the theoretical wave drag of the symmetrical 





configuration at zero lift. This drag increment was then added to all the 
theoretical results to obtain the curves of theoretical lift-drag ratio as a 
function of lift coefficient. As predicted theoretically, none of the modifi- 
cations changed the slope of the lift or moment curves. However, there were 
significant differences in the angle of attack at zero lift, pitching moment 
at zero lift, and maximum lift-drag ratio of the three configurations. As can 
be seen in figure 7, this method provides very accurate estimates of these 
quantities. 

The characteristics of these wing -body combinations are somewhat different 
from those of the isolated wing, as can be seen in figure 8. The wing referred 
to in this figure is the cambered and twisted wing of reference 7 used in con- 
figurations 2 and 3 and the wing-body is configuration 2 of figure 7- One may 
observe that the wing alone has a of approximately 0.01, whereas the body 

in combination with this wing has virtually no pitching moment at zero lift. 

In a similar manner, the lift curves of the two configurations are displaced. 
This indi cates the importance of the body on the overall properties of wing- 
body combinations and may help to explain some of the discrepancies which have 
occurred between theory and experiment on wing-body combinations designed by 
wing-alone computing programs. 


Effect of Body Size on Shape of the Minimum Drag Wing 

One of the unique features of this computing procedure is the ability to 
compute the surface of minimum drag for a given wing planform in the presence 
of a body which may be at incidence relative to the wind* To illustrate the 
importance of the body flow field a study has been made of the minimum drag 
shape of a given wing with different size bodies. The basic wing planform is 
a simple delta wing with 0 cot A = 0.3 which corresponds to a leading-edge 
sweep of 73.3° at a Mach number of ^2. The shape of the isolated wing with 
minimum drag for a fixed lift is illustrated in figure 9* Five sections are 
shown through this wing constructed with the leading edge in the x-y plane. 

The straight reference lines are the intersection of the section plane with 
the x-y plane. The same wing was then analyzed in combination with bodies of 
various sizes, one of which is shown in figure 10. The shape of the wing in 
the presence of this small body is similar to the shape of the isolated wing. 
This confirms the intuitive conclusion that a small body- should have a small 
effect on the optimum wing shape. As the ratio of body diameter to wing span 
was increased, the shape of the minimum drag camber surface for the wing-body 
combination varied considerably from the shapes shown in figures 9 and- 10* 1116 

effect on this shape of such variables as nose length and shape, incidence of 
the body relative to the free stream, and location of the wing on the body has 
not been adequately investigated at this time. For this reason, the shapes of 
these wings are not presented. 

The theoretical model upon which these calculations are based assumes 
inviscid flow. For this reason, these shapes should be viewed cautiously as 
far as actual experimental performance is concerned. 
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CONCLUDING REMARKS 


The computational method whose results are described in this paper provides 
the designer with the ability to predict the longitudinal characteristics of 
wing-body combinations of arbitrary planform, camber, twist, and thickness. By 
this procedure one may accurately estimate many important performance param- 
eters which are not available by the older techniques now in use. The technique 
is rapid and straightforward enough to be used in engineering design st ud ies of 
families of configurations. The program is written for the class of computers 
in wide use by the aviation industry today. With the optimization capability 
built into the program it should be possible to calculate low-drag wing-body 
combinations to be used as the basis for designing efficient aerodynamic 
vehicles. 
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TYPICAL PANEL LAYOUT FOR WING-BODY COMBINATION 



Figure 1 


COMPARISON BETWEEN EXACT LINEAR THEORY AND 
PROGRAM RESULTS FOR DELTA WINGS 
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Figure 2 
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COMPARISON BETWEEN EXACT LINEAR THEORY AND PROGRAM 
RESULTS FOR A SWEPT WING 
/9JR* 1.92; /9COt A *.6 




Figure 3 


LIFT-CURVE SLOPE OF DELTA WING, A =60° 
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Figure 4 













